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We analyze the effects of spatially extended periodic forcing on the dynamics of one-dimensional 
excitation waves. Entrainment of unstable primary waves has been studied numerically for different 
amplitudes and frequencies of additional sub-threshold stimuli. We determined entrainment regimes 
under which excitation blocks were transformed into consistent 1:1 responses. These responses were 
spatially homogeneous and synchronized in the entire excitable medium. Compared to primary 
pulses, pulses entrained by secondary stimulations were stable at considerably shorter periods which 
decreased at higher amplitudes and greater number of secondary stimuli. Our results suggest a 
practical methodology for stabilization of excitation in reaction-diffusion media with regions of 
reduced excitability. 



Dynamics of excitation waves in reaction-diffusion media can be altered by spatio-temporal periodic forcing. Ad- 
ditional (secondary) periodic stimulations superimposed on primary forcing may alter the primary excitation waves 
and entrain (lock) them to the period of secondary stimuli. 

Locking of primary waves to the period of secondary stimulations occurs at particular values of forcing periods 
and amplitudes. This resonant shift is characterized by Arnold tongues which determine different types of M : N 
(M > N, N > 1) locking responses as a function of amplitude and frequency of external forcing [TJI2]- 

It was found that phenomenon of locking manifests itself in different ways depending on spatio-temporal complexity 
of a particular reaction-diffusion system. For example, under the periodic forcing spatially uniform two-dimensional 
BZ reaction oscillations were transformed into standing wave type labyrinths of complex geometry pQ. It was also 
demonstrated that one-dimensional Turing patterns can be modulated by spatio-temporal forcing in the form of 
traveling wave [3]. 

Similar type of resonant behavior of Turing patterns was observed in experiments with photosensitive chemical 
reactions [1]. It was also shown that the presence of large amplitude periodic forcing in a one-dimensional bistable 
reacion-diffusion medium resulted in bifurcating of originaly stable wavefront into two counter propagating wavefronts 
[5]. In contrast, modulating the intensity and frequency of the periodic forcing controlled the trajectory and rotational 
frequency of two-dimensional spiral waves in the Oregonator model [6]. 

One of the major parameters which affect reaction-diffusion dynamics is the amplitude of a threshold stimulus (ex- 
citation threshold) required to initiate an excitation wave. While the outcome of over-threshold periodic stimulations 
has been discussed above in both experimental and theoritical settings PUS], the response to periodic sub-threshold 
perturbabtions smaller than the excitation threshold has not been analyzed. 

In this paper, we study the entrainment of marginally stable excitation waves by spatially extended periodic sub- 
threshold forcing in the Chernyak-Starobin-Cohen reaction-diffusion medium [7j. 
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Here u and v are dimensionless excitation and recovery functions, respectively, e is a small parameter and v r is the 
excitation threshold. A and ( control the rates of changes of excitation and recovery variables, respectively. 

The system of Eqs. 1 was solved numerically in a one-dimensional cable of finite length using a second order explicit 
difference scheme with zero flux boundary conditions [8]. Spatial, Ace, and temporal, At, steps used in the numerical 
integration were equal to 0.23 and 0.0072, respectively. The cable length, L, was equal to 150Ax. Primary forcing, 
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P(x,t), was a periodic function with an over-threshold amplitude Aq and period To. Primary stimuli were applied 
near the left end of the cable between x = 2Ax and x = 15Ax. As in some biological reaction-diffusion media [9] we 
used a simplified primary rate dependent threshold given by linear equation v r = a — /3T , a > 0, /3 > 0. 

Secondary forcing {Si(xi,t), i — 1,2, ...,n} was delivered at n equidistant locations between x — 40Aa; and x = 
140Aa;. The secondary stimuli were simultaneously activated after the wavefront resulted from the first primary 
stimulation arrived at the end of the cable. The amplitude and period of secondary stimuli were equal to A and T, 
respectively. 

Unless mentioned otherwise, values of parameters in all computations were e = 0.1, A = 0.4, (, = 1.2, Aq = 
1.4, a — 0.31, (3 = 0.0025 and n = 6. The duration of a pulse, Th, was measured as the time interval between 
consecutive intersections of u and v near their rest and excited states, respectiely (Fig. la). The steady state value 
of Th was computed after 80 primary stimulation periods. 

In the absence of secondary stimulation (A = 0), at long periods To the variable v has enough time to reach its 
ground state v r before the next stimulus is applied. However, as To approaches a critical limit, T en d, the system does 
not respond to every stimulus and exhibits unstable M : N (M > N) excitation blocks which occur due to incomplete 
recovery of control variable v. Figure lb compares phase portraits of the system at two values of T . At T = 60, u 
closely follows its nullcline and v almost completely recovers to v r = 0.16 as depicted by the intersection of the u-v 
nullclines. However, at To = T en( j = 30, deviations of u and v from their nullclines are quite significant, thereby v 
recovers to a value which is much higher than the corresponding threshold v r = 0.23. 

Analysis of pulse duration Th as a function of secondary frequency F = T -1 and forcing amplitude A > revealed a 
variety of entrainment regimes (Fig. 2). We found that the system did not respond to secondary stimuli at amplitudes 
which were smaller than critical values depicted by the curve with circular markers. For the amplitudes above this 
curve we observed intermediate M : M responses with M greater than one. For even greater amplitudes, above the 
upper curve with square markers, the system locked to secondary stimuli with consistent 1:1 responses. It should be 
noticed, that such locking occured in a wide range of frequencies of secondary stimulations at sub-threshold amplitudes 
which were 80% lower than amplitudes of primary stimuli. 

Spatio-temporal contours of u shown in Fig. 3a demonstrate expected unstable responses to primary stimulation at 
To < T enc i. Temporal dynamics of u and v, as well as w-contours, show 3:2 excitation blocks (Fig. 3b). However, in the 
presence of secondary stimulations such unstable responses can be entrained and stabilized by secondary driving even 
at To < T en d- Indeed, Fig. 4a demonstrates that 3:2 blocks can be transformed into stable 1:1 responses which evolve 
homogeneously in the entire cable except short segments located near the site of primary stimulation. Formation of 
these fully synchronized responses are preceded by very short (~ 0.02To) transient periods during which standing 
wave type oscillations of u rapidly saturate at constant excitation levels (Fig. 4b). 

When compared with primary forcing alone (A = 0), secondary sub-threshold stimuli facilitate development of stable 
pulses shorter than those at T en d = 30. Fig. 5a shows that secondary stimuli of higher amplitudes sustain shorter 
entrained pulses at progressively lower stimulation periods. The trend is slightly augmented for higher coefficients /3 
and lower amplitudes of secondary stimuli (Fig. 5b). 

Stabilization of the system due to secondary driving can be also achieved using greater number of secondary 
stimulation sources. Correspondingly, Fig. 6 demonstrates that due to initial two-fold increase of the number of 
secondary sources the region of stability is significantly extended towards shorter values of T en d- Further increase 
of the number of sources saturates these changes at higher rates for smaller values of the coefficient j3. Under these 
conditions excitation thresholds are lower, pulse wavefronts are narrower and, therefore, stabilization of dynamics is 
achieved at smaller values of T en d- 

In summary, we demonstrated that additional sub-threshold driving stimuli can entrain otherwise unstable primary 
reaction-diffusion waves and transform M : N excitation blocks into stable 1:1 spatially homogeneous responses 
synchronized in the entire cable. Compared to pulses resulted from primary forcing alone, pulses entrained by 
secondary stimulations were stable at considerably shorter periods. These periods decreased at higher amplitudes and 
greater number of secondary stimuli. Our results outline the possibility of entrainment of reaction-diffusion waves by 
sub-threshold additional driving and may be applied for stabilization of excitation in reaction-diffusion systems with 
zones of reduced excitability [TO] . 
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FIG. 1: Panel (a) shows excitation (solid line) and recovery (dashed line) variables at x — L/2 for the short steady state pulse. 
Left and right intersections between u and v mark the begining and the end of the pulse, respectively. Panel (b) shows phase 
portraits of the system at the same location for To = T e „d (thin-dashed contour) and To = 60 (thin-solid con- tour). Nullcinc 
of u is given by thick N-shaped line. Dashed and solid lines with intercepts at v T = 0.23 and v r — 0.16 are the v nullclines for 
To = T end and T = 60, respectively. 
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FIG. 3: Panel (a) shows the dynamics of u as a function of time, t, and spatial location, Xi = i Ax, {i = 0, 1, 150}. Panel 
(b) demonstrates temporal dynamics of u (solid line) and v (dashed line) for three consecutive cycles at To = 28, A = 0. x = — 
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(a) (b) 

FIG. 4: Panel (a) shows spatio-temporal band type synchronized patterns of u (To = T = 28) for A — 0.1 Ao. Panel (b) depicts 
temporal evolution of a spatial profile of u at three equidistant moments of time, starting from Unit = 136, St — 0.36. 




FIG. 5: Panel (a) shows dependence of Tu on T for different amplitudes of secondary stimulli, T = T. For each amplitude 

ratio, the lowest value of T corresponds to T en d- Panel (b) depicts dependence of T en£ j on — — for two values of /?. a = 0.31 

Ao 

and x = — . 
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